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Low-mass X-ray binaries in the bulge of the Milky Way 
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Abstract. We study the population of low-mass X-ray binaries (LMXBs) in the Galactic bulge using the deep 
survey of this region by the IBIS telescope aboard the INTEGRAL observatory. Thanks to the increased sensitivity 
with respect to previous surveys of this field, we succeeded to probe the luminosity function (LF) of LMXBs down 
to ~ 7 x 10 34 erg s" 1 in the 17-60 keV energy band. The slope d log N/d log L = -0.96 ± 0.20 measured in the 

10 35_ 10 37 

erg s range confirms that the LMXB LF flattens below L x ^ 10 3 erg s 1 with respect to higher 
luminosities. We discuss the origin of the observed LF flattening. We demonstrate that the spatial distribution 
of persistent LMXBs in the Galactic Center/Galactic bulge region is consistent with a model of stellar mass 
distribution that includes the nuclear stellar disk component in the innermost degree of the Galaxy. The spatial 
distribution of transient LMXBs detected in the Galactic Center region indicates an increased fraction of transient 
sources in the innermost degree of the Galaxy with respect to outer regions. 

Key words. ISM: general - Galaxies:general - Galaxies: stellar contents - X-rays:diffuse background 



1. Introduction 

Low-mass X-ray binaries (LMXBs) are binary systems 
with a compact object that accretes matter from a low 
mas s optical companion star filling its Roche lobe (see 
e.g van den Heuvel|[T975[) . Ways of formation and evolu- 
tion of LMXBs are important and not yet fully under- 
stood questions of X-ray astronomy (see discussions of 
these issues in e.g. lWebbink. Rappapqrt, fc Savoniielll983t 



iLipunov et"aill2007t iBelczvnski et al.ll2008f ) 

The shape of the luminosity function (LF) of LMXBs 
may provide important information about their lon g term 
evolution. In particular. iPostnov fc Kuranovl ( 2005 ) noted 
that the observed shape of the LMXB LF (assuming that 
the X-ray luminosity of the compact object is directly 
proportional to the mass transfer rate in the binary sys- 
tem) depends on the distribution of the masses of optical 
donor stars in LMXBs, which is in turn intimately related 
to the mechanism of removal of the angular momentum 
of the binary system. They suggested that the flattening 
of the LMXB LF at luminosit ies L x < 2 x 10 37 erg s" 1 
found in a number of studies JPrimini. Forman. fc Jones] 
19931: lGilfanovll2004t IVoss fc Gilfanovl I2006L l2007ri might 



itational braking mechanisms in LMXBs and enable one 
to select those binary systems effectively emitting grav- 
itational waves for the next generation of gravitational 
wave detectors. It is therefore interesting and important 
to check the presence of the break in the LMXB LF and 
to continue the LF to lower luminosities. 

Using LF as a tool for studying populations of X-ray 
sources in nearby galaxies is only possible with instru- 
ments that have good sensitivity and angular resolution. 
First attempts of estimating the LF of bright X-ray bina- 
ries were done with the Einstein/ H EAQ2 observatory (see 



van Spevbroeck et al. Il979t Trinchicri & Fabbiand 



19911 ) more than twenty years ago. Studying LMXB 
LF in galaxies has now become quite common due to 
the unprecedented capabilities of the Chandra X-ray 
Observatory. However, even this powerful observatory 
does not allow one to probe the faint end (L x < lfj 36-37 
erg s _1 ) of the LMXB LF in the majority of nearby galax- 
i es. This is only feasible for th e nearest galaxies, like M31 
( Trinchieri fc Fabbiand 1991 : Primini. Forman. fc Jonesl 



19931 : IVoss fc Girfanovl l2007) 



correspond to the transition from the magnetic stellar 
wind braking to the gravitational wave braking mecha- 
nism of_j^movaJ_oJ_th£_binaxy system angular momen- 
tum ( Postnov fc Kuranov| |2005) . Confirmation of this idea 
would provide additional evidence for the action of grav- 
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The Milky Way is a unique galaxy where we can de- 
tect X-ray objects with the lowest possible luminosities. 
However, because of the large angular size of the Galaxy 
and the different and often poorly known distances to 
Galactic X-ray sources, it is not easy to construct the LF 
of LMXBs in the Milky Way. First of all, one needs to 
cover a large fraction of the sky. However, focusing X-ray 
telescopes usually have small fields of view (~10 arcmin) 
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and thus cannot cover a considerable part of the Galaxy 
within a reasonable time. Even the region of maximal con- 
centration of Galactic LMXBs - the Galactic bulge - is too 
large for focusing X-ray telescopes. 

The Galactic Center/Galactic bulge region is most 
suitable for studying Galactic LMXBs because the ma- 
jority of sources are spatially concentrated and are lo- 
cated at nearly the same distance from the Sun. I Skinner 
(Il993h reconstructed the surface density map of X-ray 
sources in the Galactic bulge region using a compila- 
tion of observational campaigns by different observatories. 
The population of Galactic bulge X-ray sources was also 
studied with the exceptionally wide field-of-view ART-P 
coded-mask telescope aboard the GRA NAT observatory 



( Grebenev. Pavlinskv. fc Sunvaevl 1996[) . 

Grimm. Gilfanov. fc Sunvaevl (|2002l ) used data of the 



All- Sky Monitor (ASM) of the RXTE observatory in the 
2-10 keV energy band to construct the luminosity func- 
tion of Galacti c LMXBs with luminosities L K 10 36 
erg s . Later, iGilfanovl (|2004h studied the populations 
of LMXBs in the Milky Way and nearby galaxies and 
showed that the LMXB LF has a break at a lumi- 
nosity L x ~ 2 x 10 37 erg s , below which it signifi- 
cantly fla ttens. Similar flattening of the L MXB LF was 
found bv IPrimini. Forman. fc Joned (1993) for the pop- 
ulation of sources in the M31 galaxy. Recent studies of 
a number of galaxies performed with the Chandra ob- 
servatory confir med such a break in the LMXB LF in 
the M31 galaxy (jKong et al.ll2003t IVoss fc Gilfanovll2007t) 
and in dicated the presence o f a break for the Cen A 
galaxy (IVoss fc Gilfanov 2006). No breaks, however, were 
detected by iKim et al.1 ([2006) in the LF of sources in 
NGC3379 and NGC4278. 

In order to make a significant improvement in the sta- 
tistical study of LMXBs in our Galaxy and independently 
measure the faint end of the LMXB LF, one needs a 
large field-of-view instrument with increased (in compari- 
son with previous ones) sensitivity. It is also necessary to 
use a hard X-ray energy band to avoid contamination by 
large numbers of nearby, soft X-ray emitting stars, as was 



the c ase with the ROSAT all-sky survey (e.g. lVoges et al 



1999). This is exactly what is provided by t he coded-mask 



telesc opes of the INTEGRAL observatory ([Winkler et al 



l2003h . Deep observations of the Galaxy performed with 
INTEGRAL/IBIS over five years of its operation have al- 
lowed us to increase the sensitivity to point sources by an 
order of magnitude in comparison with the RXTE/ ASM 
survey. This allows us to probe the faint end of the LF of 
Galactic LMXBs for the first time. 

2. Sample of sources and data analysis 

To avoid problems with unknown source distances and 
consequently unknown source luminosities, we do not use 
here the whole INTEGRAL survey of the Milky Way 
and restrict ourselves to the Galactic Center/Galactic 
bulge re gion where LMXBs are strongly concentrated 
(see e.g ISkinnerl 11993k iGrebenev. Pavlinskv. fc Sunvaevl 



I l996 |: Grimm. Gilfanov. fc Sunvaevll2002T ; Lutovinov et al 
120051 ). 

The INTEGRAL observatory has spent a lot of 
time observing the Galactic Center (Galactic bulge) re- 
gion over the first five years of its operat ions (see e.g 
Revnivtsev et al.|[200l iKuulkers et al.ll2007t ). Among the 



INTEGRAL instruments, the IBIS telescope best meets 
the requirements of the present study: it has good sensi- 
tivity for detecting point sources and a wide field of view 
(~ 28° x 28°). Due to the hard energy band of IBIS (> 17 
keV), it is less sensitive to sources that have high lumi- 
nosities in the standard 2-10 keV energy band (since such 
sources have very soft spectra that do not continue to hard 
X-rays), but it is very effective in detecting low- luminosity 
objects (which typically have hard spectra). The sensitiv- 
ity achieved by now with respect to point sources in the 
Galactic Center region is ty pically 0.5 mCrab in th e 17- 
60 keV energy band (see e.g. iKrivonos et al.ll2007bl ). 

For the present study we selected all sources detected 
by INTEGRAL in the Gala ctic Center (Ga l actic b ulge) 
region listed in the catalog of IKrivonos et al. I (|2007bh and 
added INTEGRAL /IBIS data that became available since 
the publication of this catalog. These additional data al- 
lowed us to increase the accuracy of source flux measure- 
ments and to add one newly detected source, IGR J17586- 
2129, to the sample. 

We only considered those sources located within the 
elliptical region around the Galactic Center with axes 
/ < 10.7°, \b\ < 5.1° and having fluxes higher than 0.64 
mCrab in the 17-60 keV energy band, which typically 
ensures the source detection with statistical significance 
more than 8a. This conservative detection threshold pre- 
vents any problems with possible systematical uncertain- 
ties in the image reconstruction and with variations of the 
sensitivity over the region of our interest. 

Due to the anticipated concentration of LMXBs in the 
Galactic bulge, we assume that all sources that are not 
known to be high-mass X-ray binaries, cataclysmic vari- 
ables or extragalactic sources, are LMXBs. It is possible 
that a detailed study of some faint sources (typically with 
luminosities L < 10 35 erg s _1 ) from our sample will re- 
veal their non-LMXB nature. Therefore, we are probably 
somewhat overestimating the faint end of the LMXB LF. 

From the preliminary source list we have filtered out: 



1. HMXBs. High-mass X-ray binaries (HMXBs) are 
young objects, which in our Galaxy are concentrated 
to the regions of recent star formation - the spiral 
arms. The spatia l distribution of Galactic HMXBs has 
been studied by iGrimrm Gilfanov, fc Sunvaev 
in soft X-rays and by Lutovinov et al.l "i 2005 . 



(2002) 



2007) 



in hard X-rays (> 20 keV) with the INTEGRAL ob- 
servatory. Due to the low star-formation rate in the 
Galactic bulge we can anticipate that HMXBs should 
not be a dominant population of bright X-ray sources 
there. However, we might detect some HMXBs located 
in the Galactic disk, projected on the bulge in the 
sky plane. We excluded from our analysis all sources 
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that are known or supposed to be HMXBs, for exam- 
ple so-called supergiant fast X-ray transients (SFXTs) 
and X-ray pulsars. These sources include EXO 
1722-363, IGR J17407-2808, XTE J1743-363, IGR 
J17544-2619, IGR J17391-3021 (XTE J1739-302), 
AX J1749. 1-2733, AX J1749. 2-2725 and IGR 
J18027-2016. Note that most of them have a transient 
nature. 



Table 1. List of known or likely LMXB sources (persistent and tran- 
sients) in the Galactic bulge region detected by INTEGRAL on the 
time-average map. The 1 mCrab flux in the 17—60 keV energy band 
for a Crab-like spectrum corresponds to an energy flux of ~ 1.4 X 10 — 11 
erg s _1 cm~ 2 and to a luminosity 1.1 X 10 35 erg s — 1 for a distance of 8 
kpc. Fluxes in the 2—10 keV energy band reported for some sources in the 
last column were obtained from RXTE/ASM measurements or adopted 
from the literature if a reference is given. The hard X-ray fluxes of the 
sources SLX 1744-299 and SLX 1744-300, which are not resolvable with 
INTEGRAL/IBIS, were calculated from their measur ed summed flux by 
assuming their flux ratio to be 2:1 llSidoli et al.lll999fl 



# 


Name 


Flux, mCrab 


Flux, mCrab 






17-60 keV 


2-10 keV 


1 


GRS 1758-258 


59.07 ± 0.07 


25 


2 


GX 1+4 


54.14 ± 0.08 


11 


3 


GX 5-1 


48.14 ± 0.07 


972 


4 


GX 354-0 


37.57 ± 0.08 


99 


5 


IE 1740.7-2942 


31.79 ± 0.07 


10 


6 


4U 1724-30 


17.14 ± 0.07 


25 


7 


A1742-294 


13.71 ± 0.07 


18[1] 


8 


GX9+1 


16.50 ± 0.09 


532 


9 


SLX 1735-269 


11.21 ± 0.07 


14 


10 


GX3+1 


9.93 ± 0.07 


300 


11 


SLX1744-299 


5.47 ± 0.07 


9.0[1] 


12 


1E1742. 8-2853 


6.29 ± 0.07 




13 


1E1743. 1-2843 


5.20 ± 0.07 


8.1 [2] 


14 


SLX 1737-282 


3.94 ± 0.07 


5.4[3] 


15 


SLX1744-300 


2.73 ± 0.07 


5.1[1] 


16 


IGR J17254-3257 


1.72 ± 0.08 


3.5 


17 


IGR J17475-2253 


1.09 ± 0.08 




18 


IGR J 17353-3539 


0.88 ± 0.08 




19 


IGR J 17505-2644 


0.81 ± 0.07 




20 


IGR J 17585-3057 


0.81 ± 0.07 




21 


IGR J17586-2129 


0.73 ± 0.08 




22 


IGR J17448-3231 


0.70 ± 0.07 






Transients detected on the average 


i map 


1 


GRS 1741.9-2853 


3.17 ± 0.07 




2 


SAX J1747.0-2853 


3.30 ± 0.07 




3 


KS 1741-293 


4.95 ± 0.07 




4 


XTE J1739-285 


2.21 ± 0.07 




5 


GRS 1747-313 


1.53 ± 0.07 




6 


IGR J17464-3213 


22.07 ± 0.07 




7 


SLX 1746-331 


0.69 ± 0.07 




8 


IGR J 17353-3257 


1.12 ± 0.07 




9 


MXB 1730-33 


4.20 ± 0.08 




10 


IGR J17331-2406 


1.01 ± 0.08 




11 


XTE J1720-318 


1.73 ± 0.08 




12 


IGR J17597-2201 


4.62 ± 0.08 




13 


XTE J1817-330 


7.01 ± 0.09 




14 


4U 1746-37 


2.77 ± 0.09 




15 


4U 1705-32 


2.30 ± 0.09 




16 


IGR J17091-3624 


4.75 ± 0.10 




[1] 


- Sidoli et all ((1999), 


2] - Porquet et al. 2003, [3] - in't 



Zand et al. 2002 



2. Extragalactic. We also excluded extragalactic 
sources. Within the region of our study there are only 
two known extragalactic sources - active galactic nu- 
clei GRS 1734-292 and IGR J17488-3253. According 
to the all-sky averaged AGN numb er-flux function in 
hard X-rays (jKrivonos et al.ll2007bl ) we should detect 
~ 2 AGNs above the adopted flux limit of 0.64 mCrab 
within the area of our study (~ 171 sq. deg). This in- 
dicates that most likely we already know of all the ex- 
tragalactic sources with fluxes of > 0.64 mCrab within 
the considered field. The probability that there are ad- 
ditional 2 or more unidentified extragalactic sources in 
this region is ~0.14 and we do not expect unknown 
extragalactic sources to considerably contribute to the 
measured number-flux function of sources within the 
bulge region. 

3. Other sources We removed cataclysmic variables, de- 
tected in the considered region. We removed from our 
sample the source IGR J17402— 3656, which is asso- 
ciated with the nearby open star cluster NGC 6400 
and may be composed of multiple point sources; the 
source IGR J 17456— 2901, which is similarly associated 
wit h the collective emission of the nuclear stellar clus- 
ter ( Belanger et al. 20061 : Krivonos et al.1 [2007a ); and 
the source IGR J17475— 2822 repre senting the emis- 
sion o f the molecular cloud Sgr B2 (jRevnivtsev et al 
2004bh . 



2.1. Transients 

A significant fraction of Galactic LMXBs are transients. 
During outbursts their luminosities can change by more 
than 6-7 orders of magnitude (e.g. iTanaka fc Shibazakil 
1996). The "on" state of bright X-ray transients typically 
lasts several weeks, but due to their extreme brightness at 
that time they may still be detected on a sky map averaged 
over several years even if they were undetectable most of 
the time. 

It is usually considered that the observed X-ray bright- 
ness variations of LMXB transients are caused by a devel- 
oped mass transfer instabil i ty in the accretion disk (e.g. 
Mever fc Mever-Hofmeisterl 119841 lLasotal l200lh . In this 
case, the X-ray luminosity observed from such a source 
at any given time does not provide us information about 
the global mass transfer rate in the binary system, which 
is of interest for us. The global mass transfer rate could be 
determined by averaging the X-ray luminosity of a tran- 
sient source over a period that is much longer than the 
accretion disk mass accumu lation time; however, this time 
scale might be too long (e.g. Chen. Shrader. fc Livio|[l997 ) 
to be testable by available X-ray observations. The inclu- 
sion of transient sources with their fluxes averaged over 
some arbitrary period (not related to any physically mo- 
tivated time scale) might distort the true distribution of 
mass transfer rates, in which we are interested. 

Therefore, in our analysis we treated the transient 
sources in two ways: 1) we included all sources detected at 
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Fig. 1. Map of the Galactic Center region in the 17-60 keV 
energy band obtained by INTEGRAL/IBIS. The ellipse 
encloses the region of our study. Sources used in the anal- 
ysis of the LMXB LF are marked by circles. Contours are 
isophotes of the 4.9/ito surface brightness of the Galaxy 
(COBE/DIRBE) demonstrating the bulge/disk structure 
of the inner Galaxy. 



the 5-year averaged INTEGRAL/IBIS map of the Galactic 
Center, and 2) we filtered out all transient sources from 
our sample. In the latter case, we adopted the following- 
selection criterium for separating transient sources from 
persistent ones: a transient should have the ratio of maxi- 
mal 3-day averaged flux to the all-time mean flux of more 
than 5.0. We note that transient nature is difficult to as- 
certain for weaker sources because they become detectable 
only after accumulating a large part of the total exposure 
time. Therefore, our resulting sample of transient sources 
is likely significantly incomplete at low fluxes. 

The final list of accepted sources is presented in Table 
[TJ The map of the Galactic bulge region with the sources 
from Table Q] is presented in Fig. [TJ 

2.2. Distances to the sources 

The distances to the majority of our sources are unknown. 
However, assuming that LMXBs are distributed in the 
Galactic bulge similarly to ordinary stars, we can calculate 
the density distribution of LMXBs in the bulge and pre- 
dict the spread of values of the flux to luminosity conver- 
sion coefficient for our sources. This determines the typical 
systematic error we are making in estimating source lumi- 
nosities by assuming that all sources are at the same, 8 
kpc, distance from the Sun. 

To this end, we adopted a mass model of the Galaxy 
consisting of the Ga lactic bulge, Galactic disk and nu - 



clear stellar disk (e.g. Launhardt. Zvlka. fc Mezgerll2002l ). 
We did not consider the nuclear stellar cluster, which 
is important in the innermost ~ 10' of the Galaxy, be- 
cause its size is smaller than the angular resolution of 
INTEGRAL/IBIS. 

The stellar density in th e bulge was assumed to follow 
model G3 from lDwek et al. ( 1995). We adopted the mass 
of the bulge to be 1.3 x 10 10 M g (|Dwek et al.lfl995h . 




-0.5 0.5 

Logarithm of lum. correction factor (log-joCF) 

Fig. 2. Probability density of the luminosity correc- 
tion factor (the probability density is normalized so 
that / P(log 10 CF)d(log 10 CF) = 1). Eighty per cent of 
sources in the studied region of the Galaxy are expected 
to fall into the interval shown in gray. 



The stellar density in the Galactic disk: 



Pdisk(r,z) = po.diskexp 



?"disk 



Zdisk 



wh ere ryiisk = 2.2 kpc, r m = 2.5 kpc, Zdisk = 130 
pc (jRevnivtsev fc Sazonovl 120071 ) . The total mass of the 
Galactic stellar disk was taken to be 2.5 x 10 10 Mfr 



i . 



The nuclear stellar disk (NSD) was assumed to have 
the density distribution (r and z are measured in parsecs) 



PNSD = /W- Q e-l z l/ 2d , 



where Zd = 45 pc. At r < 120 pc the slope a = 0.1 
and the normalization constant pd = 300M Q /pc 3 . At 
120 pc < r < 220 pc, a = 3.5, and at r > 220 
pc, a = 10 (jLaunhardt. Zvlka. fc Mezger 2002 ). and we 



adjust the constant pd so that the density distribution 
is continuous. The total adopted mass of the NSD is 
1.4 x 10 9 M(T) . In reality this quantity is u ncertain by some 
50% (|Launhardt. Zvlka. fc Mezgerll2002t) . 

Based on this model, we determined the distribution 
of stars over the distance from the Sun in the region of our 
study (ellipse \l\ < 10.7°, |6| < 5.1°). The ratio (D/8 kpc) 2 
defines the luminosity correction factor CF. The proba- 
bility density of this correction factor can be calculated 
from the probability density P(D) of the source distance, 
which can be found by integrating the distribution of stel- 
lar mass p*(l, b, D) in a given direction (I, b) over the solid 
angle £1 of our study: 



P(D) 



f nP *(i,b,D)D 2 dn 

LJ"d=T P*(l,b,D)D*dQdD' 
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Fig. 3. Typical broadband spectra of LMXBs. The right 
axis shows the logarithm of the luminosity, assuming a 
distance of 8 kpc to all sources. 

The probability density of CF is presented in Fig. [2] 
It turns out that for 80% of sources log 10 CF < 0.3. 
Therefore, assuming that all bulge sources are located at 
the same distance (8 kpc) from the Sun we can predict 
the luminosities of most sources to within a factor 2, this 
uncertainty being comparable to expected variations in 
source fluxes. 

2.3. Spectral corrections 

Previously, for constructing the luminosity function of 
LMXBs in galaxies one usually used instruments operat- 
ing in the s tandard X-ray energy band, 2-10 keV, or sim- 
ilar b ands (jGrimm. Gilfanov. &: Sunvaevl 120021 : iGilfanov 
2004) . The INTEGRAL/IBIS sky survey is done in the 
hard X-ray band, 17-60 keV. Therefore, to compare the 
LMXB LF derived here with those reported before we 
should take into account the spectral shapes of LMXBs 
in the broad energy range 1-100 keV. 

It is well-known that the brightest LMXBs emit 
most of their luminosity in the standard X-ray band 
(~2-10 keV), while less luminous objects emit similar 
amounts of energy in the standard and hard X-ray (~10— 
100 keV) bands. Bright LMXBs usually have spectra 
with an exponential cu t off at energ i es of ~6-10 keV (e.g 



Rappaport et al.| 1969; 



Toor et al.l 



gies ot r> 



Mitsuda et al 



19841 IPavlinskv. Grebenev. &; Sunvaevl Il994h and essen- 
tially do not contin ue into the hard X-ray band (e.g. 
Gilfanov et a"Dll993l) . In contrast, dim LMXBs are char- 



acterized by a power-law spectral shape between 2 



and 60 keV ([Mitsuda et al. 1989; Grebenev et al 



Barret et all 1 19911: IZhang et al.Hl996t iBarret et all 



1991 



2000) 



We can therefore expect that the LMXB LF constructed 
in the 17-60 keV energy band should strongly differ from 
that derived in the standard X-ray band. 

To demonstrate this, we constructed the broadband 
spectra of Galactic bulge sources with different luminosi- 
ties using all available instruments - INTEGRAL/IBIS, 
RXTE/PCA, Swift/XRT, ASCA/GIS. This allowed us to 
cover the wide energy range 1-100 keV (see Fig. [3]). The 
LMXB GS 1826-24 formally lies outside the elliptical re- 
gion of our study, but since its distance from the Sun 
is sim ilar to those of bulge sources (e.g 



in 't Zand et al 



1999) and its broadband spectrum has very high statistics 
this spectrum is also shown in the figure. 

One can see dramatic changes in the broadband be- 
havior of source spectra at a threshold luminosity of 
~ few x 10 37 erg s _1 . At high luminosities, most of the 
energy is emitted below 17 keV. At low luminosities, the 
spectral shape in the 2-60 keV band is close to a power 
law with photon index L ~ 1.9-2 and similar fluxes are 
emitted at 2-10 keV and 17-60 keV. 

The time-averaged fluxes of most of the brightest 
sources in the 2-10 keV band are directly measured by 
the All-Sky Monitor aboard RXTE and we present these 
measurements in Table [1] Based on these data and Fig. [3J 
we adopted the following spectral correction factors for 
the subsequent analysis: 



L 2 _ 



2-10 koV 



Li 



7-60 keV 



40 if L 2 _io kcv > 2 x 10 37 erg s" 1 
1.25 if L 2 -io kcv < 2 x 10 37 erg s" 1 



In comparing below the LMXB LFs in the 2-10 keV 
and 17-60 keV energy bands, we mostly used direct 
RXTE/ ASM measurements to determine the 2-10 keV lu- 
minosities for sources with £2-10 kcv ^ 10 36 erg s _1 , and 
assumed that those sources with luminosities below this 
RXTE/ ASM completeness limit are in the hard spectral 
state. 

It is important to note that the spectral correction 
factors introduced above may not be applicable to rare, 
extremely luminous LMXBs such as the Galactic source 
GRS 1915+105, which is more powerful that any source 
in our Galactic bulge sample and is apparently in a highly 
unusual spectral state. 

3. Results 

3.1. Luminosity function of faint LMXBs 

The luminosity functions of all and only persistent LMXBs 
detected by INTEGRAL/IBIS in the Galactic bulge region 
are shown in Fig. 2] These LFs are defined so as to give 
the number of sources within the studied region containing 
1.65 x 10 10 M Q of stars. 

It can be seen that the LF of weak LMXBs 
(ii7_ 60 kcv < 10 36 - 5 erg s^ 1 ) is relatively flat. The 
LF of all LMXBs can be approximated by the function 
L 35 dN/dL 35 = (12±4)L~ (al3±0 - 13) in the range 10 34 7 
erg s" 1 < L17-60 kcv < 10 36 - 5 erg s" 1 , and the LF of 
peristent sources can be approximated by the function 
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35 36 37 

lo 9 l 17-60 kevAerg/sec) 

Fig. 4. Luminosity function of LMXBs detected by 
INTEGRAL/IBIS in the Galactic bulge region. The LF 
of persistent sources is shown by solid crosses, and that 
of all sources (including transients) by dashed crosses. 
The thick gray cross represents an estimate of the num- 
ber of LMXBs with extremely high luminosities based 
on the one such source in the Galaxy, GRS 1915+105 
(Li7_6o kcV ~ 5 x 10 37 erg s _1 ). This estimate takes into 
account that the bulge region covered by the present study 
contains only 16% of the total stellar mass of the Galaxy. 
The thic k gray lin e is th e analytic model of the LMXB 
LF from IGilfanovl \2004) normalized to 1.65 x 1O 1O M 
- the stellar mass within the considered region, recalcu- 
lated into the 17-60 keV energy band using the introduced 
spectral correction factors and multiplied by an additional 
factor of 0.4. The thick blue line is the same LF model 
convolved with the distribution of luminosity correction 
factor, which is important for sources in our sample due 
to unknown values of exact distances (see Fig[2|) 

L 35 dN/dL 35 = (4.7 ± 2.2)45° 4±0 ' 2) - Here, L 35 ~ is the 
hard X-ray (17-60 keV) luminosity in units of 10 35 erg s _1 . 

There are no objects with Lnso k c v > 10 37 erg s _1 
in our sample. Consequently, there is an abrupt drop in 
the number density of very bright sources as can be seen 
in Fig. |4] In fact, there is only one non-transient (but 
strongly variable) source in the whole Milky Way that has 
a higher luminosity: GRS 1915+105 with Li7_60 keV ~ 
5 x 10 37 erg s _1 (this value was calculated from the time 
averaged flux of the s ource of 261 mCrab in the energy 
band 17-60 keV from iKrivonos et~al . 2007b, assu ming a 

Since 



Harlaftis fc Greine 



r 2004). 



11 kpc source distance, 
such ultraluminous sources can be easily detected with 
INTEGRAL throughout the Galaxy, we can estimate their 
expected number in the Galactic bulge region studied here 
by taking into account that this region contains 16% of the 



in 

ro 



in 
ro 




0.1 r 



36 37 38 

'°9 L2-10 kevAerg/sec) 

Fig. 5. Luminosity function of persistent sources in the 
2-10 keV energy band. The solid line s hows the analytic 
form of the LMXB LF from IGilfanovl (|2004l) rescaled to 
1.65 x 10 10 M Q - the stellar mass within the region stud- 
ied here, and multiplied by 0.4. The shaded ar ea repre- 
sents the spread of LF between different galaxies (|Gilfanov 
2004h 



total stellar mass of the Galaxy. This estimate is shown in 
Fig. 0]by the thick gray cross. 

In order to compare the 17-60 keV LF derived here 
with the a nalytic approxim ation of the 2-10 keV LF pre- 
sented by IGilfanovl (|2004r ) , we recalculated the last one 
into the 17-60 keV energy band using the spectral correc- 
tion factors described above. This predicted LF, rescaled 
to the stellar mass contained in the considered region 
(1.65 x 10 10 M Q ), is shown in Fig.[4]by the thick gray line. 
It was additionally multiplied by a factor of 0.4 in or- 
der to match the normalization of the LF measured with 
INTEGRAL in the Galactic bulge. We note that a similar 
indication of underabundance of LMXBs in the Milky Way 
as compared to other galaxies was previousely reported by 
Gilfanovl (|2004 . 

As is apparent from Fig. 2J transient sources provide a 
significant contribution to the LMXB LF. If we describe 
the LF by the function L^^dN/dL^ = const in the lu- 
minosity range 10 34 7 -10 36 5 erg s _1 , the normalization 
increases by a factor of ~ 1.7 after adding the transient 
sources. 

In Fig. [5] we present the luminosity function of per- 
sistent LMXBs in the standard X-ray band (2-10 keV), 
which we constructed using direct RXTE/ASM source 
flux measurements for brighter sources (see Table [lj and 
our assumed spectral correction factors for weaker sources. 

We should note that, despite the low star- formation 
rate in the Galactic bulge, some bulge HMXBs may be 
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present in our sample. This can be especially important at 
the lowest luminosities because the LF of HMXBs is much 
stee per than that of LMXBs below ~ 2 x 10 37 erg s~ 1 (see 



e.g. iGrimm. Gilfanov. fc Sunvaevl 12002 : iGilfanovl 120041 1. 
For example, even for a star formation rate as low as 
M* = 0.03M Q /year, the number of faint HMXBs will be 
comparable to that of all similarly faint sources in the 
region of our study. Therefore, the true number of faint 
LMXBs will be even lower than suggested by the derived 
luminosity function. 

Summarizing all of the above, we can conclude that 
irrespective of the inclusion of transient sources the shape 
of the LF of Galactic bulge LMXBs is broadly compat- 
ible with the L MXB LF averaged over nearby galaxies 
(|Gilfanovil2004[ ). which exhibits a significant flattening at 
luminosities L2-10 keV < 2 x 10 37 erg s _1 . 



d log N/d log M = —a and L x oc M 13 , then the resulting 
slope of the LF will be 7 = (a + [3 — l)/(3. Assuming 
that the true mass transfer rate distribution has the slope 
a = 1.8-2.0, a s seen at luminosities log£ x ^37.5-38.5 
(IGilfanovl libol . and that = 3.3 (jBeckert fc Duschl 
20021) . the slope of the LF will be 7 ~ 1.2 - 1.3. Such 



a flattening of the LF is significant but not strong enough 
to explain the observed behavior of the LF. 

The majority of the considered X-ray sources are likely 
accreting neutron stars. They have solid surfaces and the 
simple advection dominated model, in which the accreted 
mass flows under the black hole horizon without creation 
of the X-ray emission, is not directly applicable to them. 
However, some other mechanisms of reducing the X-ray 
emissivi ty of the accretion system might be consi dered 
(see e.g. lTorbett|[l983 iBlandford fe Begelmanlll999f) . 



3.2. Possible interpretations of the break in the LF 

The flattening of the LMXB LF at luminosities L x ~ 10 37 
erg s _1 could be due to various reasons: i) a break in the 
distribution over mass transfer rates from the companion 
star, ii) a change of the correspondence between the mass 
transfer rate in the binary system and the mass accretion 
rate onto the compact object, possibly associated with the 
spectral state transition. We now briefly comment on each 
of these possibilities. 



3.2.1. State transition 

It is somehow surprising that the break in the LMXB LF 
occurs almost exactly at the (broadband) luminosity at 
which the sources experience a spectral state transition. 
At high luminosities, practically all sources are in the so- 
called soft /high spectral state and emit most of their lumi- 
nosity in the optically thick regime. At low luminosities, 
the sources are in the so-called low/hard spectral state, 
emit in the optically thin regime and have hard spectra. 

It is usually assumed that in the soft/high spec- 
tral state the X-ray luminosity is directly proportional 
to the mass accretion rat e in the binary system (e.g. 
Shakura fc Svunvaevlll973l ). The dependence of X-ray lu- 
minosity on mass accretion rate in the low/hard state 
is not clear. In the frame work of the radiatively ineffi- 
cient accretion flow models ( Ichimaru 1977t|ToYbetd[l984 : 



3.2.2. Change of the binary braking mechanism 



Naravan fc Yll994lBlandford fc Begelmanlll999l ). the ac 
cretion flow X-ray luminosity in the hard state may be 
a non-linear function of the mass accretion rate at the 
outer edge of the accretion flow: £ x oc M' 3 , where (3 > 1. 
Some observatio nal data favor such a dependence (e.g. 



Gallo et al.ll200l 



In this case, even if the distribution of mass transfer 
rates in binary systems has no peculiarities at the rates 
corresponding to X-ray luminosities ~ few x 10 37 erg s _1 , 
one would still observe a break in the X-ray luminosity 
function due to the change in the L x — M relation. If the 
underlying distribution of mass transfer rates has the slope 



Following iPostnov fc Kuranovl (|2005f ). the break in the 
LMXB LF can be interpreted as the luminosity at which 
the mass transfer rate in the binary systems starts to 
driven by the gravitational braking mechanism rather 
than by magnetic stellar wind braking. 

If this is indeed the case, we can anticipate that the 
majority of such faint LMXBs should have very small 
mas ses of optical companions (M opt £5 0.2 — 0.4M Q , 
see IPostnov fc Kuranovl [2005) and short orbital peri- 
ods (< 2 — 3 hours), and be very faint in the op- 
tical and near- infrared spectral bands: the observable 
brightness of sources could be as low as itir > 20 
even without any interstellar extinction (see estimates 
of the optical brightness of com pact binary systems in 
van Paradiis fc McCrintocklll994l ). These predictions can 
be directly checked when accurate localizations of such 
X-ray sources will be achieved with the help of X-ray tele- 
scopes. Positional accuracy of the order of few arcseconds 
or less is necessary in this case, because the spatial density 
of such faint optical objects (jur > 20 — 22) is very high 
in the direction of the Galactic bulge. 

3.3. Distribution of LMXBs in the Galactic Center 
region 

3.3.1. Persistent sources 

An extensive study of LMXB distributions in the 
Milky Way and in other galaxies showed that 
their number density cl osely traces the stellar 



mass d istrib ut ion (e.g. iGurskv fc Schreierl 
Skinnerl 1993 : Grebenev. Pavlinskv. fc Sunvaev 



1975 



1996 



Grimm. Gilfanov. fc Sunvaevl 120021 : IGilfanovl \2004 ). A 
rough comparison of the surface density of LMXBs 
in the Galactic bulge r egion w it h th at of stars 
was previously done by I Skinner (1993) based on 
arc hival observations with differen t instr uments and 
by iGrebenev. Pavlinskv. fc Sunvaev ( 19961 ) based on 
a smaller sample of sources detected with the ART-P 
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Distance from GC, deg Distance from Sgr A* (longitute) 



Fig. 6. Cumulative number of persistent LMXBs in the 
Milky Way bulge seen by INTEGRAL/IBIS within the 
circular region of a given projected radii. The solid curve 
shows the number of sources predicted by the Gal axy mass 
model described above and the LMXB LF of Gilfanov 
(|2004 ) (multiplied by 0.4). The contributions of the 
Galactic bulge and nuclear stellar disk components are in- 
dicated by the short- and long-dashed lines, respectively. 
For the model curve we assumed the luminosity limit of 
10 35 erg s _1 . Positions of sources in the central region of 
the bulge are labeled. 

telescope aboard the GRAN AT observatory. Now, using 
the sample of sources detected by INTEGRAL over five 
years of its operation we can make a somewhat better 
comparison. 

The cumulative angular distribution of the number 
of persistent LMXBs detected by INTEGRAL in the 
Galactic bulge region as a function of projected angular 
distance from the Galactic Center is shown in Fig. [5] by 
the histogramm. The solid curve in this figure shows the 
cumulative stellar mass within given projected radii. To 
match the number of observed sources at distances < 5.1° 
from the Galactic Center with the galactic mass model 
described above (the total stellar mass contained in this 
region is » 1.2 x 10 10 M Q ) and the average LMXB LF of 



Gilfanov (|2004h . the latter was multiplied by a factor of 
0.4. 

Interestingly, the radial distribution of sources shows 
an indication on an excess of the number of LMXBs in- 
side the innermost ~ 2° with respect to that expected 
from the distribution of stars in the Galactic bul ge (note 
that a similar effect was previousely observed by I Skinner 
19931) . The excess itself has low statistical significance - 
instead of the predicted ~ 1 source within the innermost 
0.6° around Sgr A* we see three sources, but it is consis- 



Fig. 7. Positions of all sources detected by IBIS/ISGRI 
within the inner 1.3° around the center of the Galaxy 
(Sgr A*). Dashed (black) curves are contours of the iso 
surface density of the stellar mass in the model of the 
Galaxy, which includes the bulge, the nuclear stellar disk 
and nuclear stellar cluster (which is important only within 
^10 ' around Sgr A*). Dotted curves are iso surface den- 
sity contours of the nuclear stellar disk component only. 
Filled (gray) circles show positions of persistent sources, 
open circles - positions of transient ones. 

tent with the presenc e of LMXBs from the so-called nu- 
clear stellar disk (e.g. Launhardt. Zvlka. fc Mezgerl 120021 ) 
(its contribution is shown by the long-dashed line). This 
allows us to tentatively suggest that at least the innermost 
sources, namely 1E1742.9-2853 and 1E1743-2843, reside 
in the nuclear stellar disk. 

In order to visualize the distribution of sources in the 
innermost region of the Galaxy (see Fig. [7|), we present 
the positions of detected LMXBs within 1.3° around Sgr 
A* with overlayed contours of iso surface density of stars 
in the region for the different components of the Galactic 
stellar population. 

It is necessary to note that the best-fit LF normaliza- 
tion for the Galactic bulge population of LMXB s is ~ .4 



of th e average over a number of nearby galaxies ( Gilfanov 
20041 ). Nevertheless, the Galactic bulge normalization is 
within t he spread of be st-fit values obtained for different 
galaxies (|Girfanovll2004[ ). 



3.3.2. Transients 

Using the large time span of the INTEGRAL observations 
we have the ability to effectively determine the "transicnt- 
ness" of sources and study these transient sources sepa- 
rately. In particular, it is interesting to analyze their spa- 



Revnivtsev M. et al.: LMXBs in the bulge of the Milky Way 



9 



10 

CD 



4— 

o 

1— 
CD 
-D 

E 



■ 










297 


1 1 

o 


■ 

- 










2802 
59-285 


\ 1736- 


EJ1751 

1 












7419- 


cr 












U~l 
CO 

r\i 

6 
r< 

S KS 1741-293 


—> X 
EC 

- r 






■ 
























aj / 










GRS 1741. 


5-2853 












AX J1 745.6-2901 
















.... 1 / 


















Distance from GC, deg 

Fig. 8. Cumulative number of transient LMXB sources in 
the Galactic Center region, as seen by INTEGRAL/IBIS. 
The solid line represents the adopted mass model of the 
Galaxy normalized to the number of sources within 5.1° 
of the Galactic Center. The positions of sources in the 
innermost region of the bulge are labeled. 



Table 2. List of transients seen by INTEGRAL/IBIS over 
the period 2003-2007 



# 


Dist from Sgr A*, deg 


Name 


1 


0.02 


AX J1745. 6-2901 


2 


0.17 


GRS 1741.9-2853 


3 


0.33 


SAX J1747.0-2853 


4 


0.36 


KS1 741-293 


5 


1.27 


IGR J17419-2802 


6 


1.36 


XTE J1739-285 


7 


1.51 


GRS 1736-297 


8 


2.02 


XTE J1751-305 


9 


2.54 


GRS 1747-313 


10 


3.22 


IGR J17464-3213 


11 


4.29 


SLX 1746-331 


12 


4.48 


IGR J 17353-3257 


13 


4.67 


XTE J1807-294 


14 


5.11 


MXB 1730-33 


15 


5.64 


IGR J17331-2406 


16 


6.18 


XTE J1720-318 


17 


7.14 


A 1744-361 


18 


7.68 


IGR J17597-2201 


19 


7.95 


XTE J1817-330 


20 


8.10 


4U 1746-37 


21 


8.56 


XTE J1818-245 


22 


8.57 


4U 1705-32 


23 


8.68 


4U 1745-203 


24 


10.60 


IGR J 17098-3628 


25 


10.66 


IGR J17091-3624 



tial distribution in view of the predictions of some theories 
that the densest central regions of galactic bulges might 
contain increased numbers of tra nsients due to their en - 
hanced dynamical formation (e.g. IVoss fe Gilfanovll2007l ). 

The list of transient sources detected by INTEGRAL 
within the region of our study is presented in Table [2] 
along with their distances from the Galactic Center. The 
cumulative number of these sources as a function of their 
distance from the Galactic Center is shown in Fig. [H The 
cumulative mass of stars within the considered region, nor- 
malized to the number of transients at < 5.1° is shown by 
a solid line. One can note that there is a slight indica- 
tion that transient sources are more concentrated towards 
the Galactic Center than the persistent ones, but the sta- 
tistical significance of this effect is low. We shoul d note 
that a similar effect was observed bv lSkinnerl <| 19931) . since 
the majority of the sources considered in that work were 
transients detected with different instruments over a long 
period. If the effect of the increased surface density of 
transient sources in the center of the Galaxy is real, this 
might be related to the mechanism of dymanical formation 
of LMXBs proposed bv lVoss fe Gilfanovl (|2007h . 



4. Conclusions 

We analyzed the luminosity function and spatial distribu- 
tion of low-mass X-ray binaries in the central part of our 
Galaxy (Galactic bulge) using data from the INTEGRAL 
observatory in the 17-60 keV energy band and found that: 



No sources in the considered region demonstrate 
hard X-ray (17-60 keV) luminosities higher than 
£17-60 keV ~ 10 37 erg s" 1 . This happens because: a) 
sources with high bolometric luminosities have soft 
spectra with only a small portion of the bolometric 
luminosity emitted in the hard X-ray band, b) those 
sources with hard spectra do not have high bolomet- 
ric luminosities. In fact, in our Galaxy there is only 
one source with hard X-ray luminosity as high as 
£i7-60 kcv ~ 5 x 10 37 erg s" 1 - GRS 1915+105, which 
is pecular in many respects. 

The luminosity function of persistent LMXBs is fairly 
flat at the faint end (Li7_6o 10 37 erg s _1 ): the slope 
is dlogA/dlogi = -0.96 ±0.20. The shape of the LF 
of the Galactic bulge LMXBs in the standard X-ray 
band is consistent with that of t he LMXB LF aver- 
aged over several nearby galaxies (|Gilfanovll20"ol and 
thus confirms the p resence of the break at ~ 10 37 
pro posed ear li er by Primini. Forman. fc Jonea (|l993 ) 
and lGilfanovl (|2004 . 

The origin of the observed break in the LF is not 
clear. One of the possibilities is that it is connected 
with a change of the mechanism of removal of the 
angular momentum of the binary system: from mag- 
netic stellar wind braking to graviational braking (see 



Postnov fc Kuranovll2005l ) . We propose that some ad- 



ditional flattening of the LMXB LF might be provided 
by the change of the mass accretion rate-X-ray lumi- 
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nosity relation expected in radiatively inefficient accre- 
tion ffow models. 

4. The cumulative angular distribution of the persistent 
LMXBs in the Galactic center region traces the cumu- 
lative number of stars in the model consisting of the 
Galactic bulge and nuclear stellar disk components. 
The data favor (though with low statistical signifi- 
cance) the presence of the nuclear stellar disk com- 
ponent in the distribution of LMXBs. 

5. The cumulative angular distribution of the transient 
LMXBs slightly differs from that of the persistent 
sources, indicating some increase of the fraction of 
transient sources in the innermost regions of the 
Galaxy. However, the statistical significance of this in- 
crease is low. 
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